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ABSTRACT 

The Fermi Gamma- Ray Space Telescope reveals a diffuse inverse Compton signal in the inner Galaxy 
with a similar spatial morphology to the microwave haze observed by WMAP, supporting the syn- 
chrotron interpretation of the microwave signal. Using spatial templates, we regress out tt*^ gammas, 
as well as IC and bremsstrahlung components associated with known soft-synchrotron counterparts. 
We find a significant gamma-ray excess towards the Galactic center with a spectrum that is signif- 
icantly harder than other sky components and is most consistent with IC from a hard population 
of electrons. The morphology and spectrum are consistent with it being the IC counterpart to the 
electrons which generate the microwave haze seen at WMAP frequencies. In addition, the implied 
electron spectrum is hard; electrons accelerated in supernova shocks in the disk which then diffuse a 
few kpc to the haze region would have a softer spectrum. We describe the full sky Fermi maps used 
in this analysis and make them available for download. 
Subject headings: diffuse radiation — microwaves — gamma-rays 



1. INTRODUCTION 

The most detailed and sensitive maps of diffuse mi- 
crowave emission in our Galaxy have been produced by 
the Wilkinson Microwave Anisotropy Probe (WMAP). 
An analysis of the different emission mechanisms in these 
maps uncovered a microwave "haze" towards the Galac- 
tic center (GC) tha t has roughly sph erical morphology 
and radius ^ 4 kpc (iFinkbein ei* '200 4aD. Since its discov - 
erv. lFinkbeine^ (|2004b[ ) andiDoblcr fc Finkbeineil (j2008l ) 
have argued that the microwave haze is hard synchrotron 
emission due to the fact that alternative hypotheses such 
as free- free (thermal bremsstrahlung of the ionized gas) , 
spinning dust, or thermal dust have difficulty explaining 
its morphology, spectrum, or both. 

However, the 23-33 GHz spectrum of the haze syn- 
chrotron is harder than that expected from diffused elec- 
trons originally accelerated by su pernova shocks in the 
plane (jPobler fc Finkbeineij|200"8[) . Some variation is ex- 
pected in the synchrotron spectrum, but generally in the 
sense that it should be softer at higher frequencies since 
the electrons lose energy preferentially at high energies as 
they diffuse from their source. While there are significant 
uncert ainties in the spectrum of th e haze (see the discus- 
sion in lDobler fc Finkbeineil[2008[ ) . the data require that 
the diffused spectrum be roughly as hard as the expected 
injection spectrum from first-order Fermi acceleration at 
supernova shock fronts (number density dN/dE oc E~^). 
That is, if the electrons were produced in shocks in the 
disk, then they would have to have undergone no diffu- 
sive energy losses over a ^ 47r/3 (4 kpc)"^ volume, which 
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seems unlikely. Furthermore, there is significant emission 
in WMAP 23, 33, and 41 GHz bands from electrons that 
were generated in SN shocks; this emission has a very 
disk-like morphology (and softer spectral index which is 
consistent with shock acceleration), while the haze has a 
more spherical morphology.^ 

Together, the haze spectrum and morphology imply ei- 
ther (1) a new class of objects distributed in the Galactic 
bulge and largely missing from the disk; (2) significant 
acceleration from shocks several kpc off the plane towards 
the Galactic center; or, perhaps most intriguingly, (3) a 
new electron component from a new physical mechanism. 
The claim that novel physics or astrophysics is required 
to explain the WMAP data is called the haze hypothe- 
sis to distinguish it from the two null hypotheses: (1) 
that the microwave haze is not synchrotron, but rather 
some combination of free-free and spinning dust; and (2) 
that the haze is synchrotron, but the electron spectrum 
required is not unusual. In this work we do not address 
the origin of the electrons, but instead consider what the 
data from the Fermi Gamma-ray Space Telescope"^ imply 
for their existence. 

Electron cosmic rays at 10—100 GcV primarily 
lose energy in the diffuse interstellar medium by pro- 
ducing synchrotron microwaves and inverse-Compton 
(IC) scattered gammas. Synchrotron losses are pro- 
portional to magnetic field energy density, Ub = 
B'^/Stt, while IC losses are proportional to the inter- 
stellar radiation field energy density, Uisbf, in the 
Thomson limit, and less in the Klein-Nishina limit. 
Bremsstrahlung off the ambient gas also occurs, but is 
expected to be sub-dominant in the regions of inter- 
est. Therefore, the best test of the haze hypothesis 
is to search for IC gammas in the Fermi data, which 
was studied in the context of dark matter signals by 

® Hereafter "spherical" is taken to mean "not disk-like" - if any- 
thing, the ha ze is non-spherical in the direction perpendicular to 
the disk (see jT2l l. 

^ See http://fermi.gsfc.iiasa.gov/ssc/data/ 
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Fig. 1. — Top and bottom left: Full sky Fermi 7-ray maps in various energy bins. The mask includes the 3-month Fermi point source 
catalog as well as the LMC, SMC, Orion-Barnard's Loop, and Cen A. Bottom right: Exposure time map with our mask overlaid and 
stretched to — 100% of the peak exposure time. The variation in the exposure is a small modulation - even setting it to unity does not 
change our qualitative results. 

strain hypotheses about their origin.^ 

In §2 we briefly review the Fermi data, describe our 
map-making procedure, and display full-sky maps at var- 
ious energy ranges. In §3 these maps are analyzed with 
template correlation techniques and resultant residual 
maps and spectra are shown. Finally, §4 presents our in- 
terpretation of the signals and discusses potential sources 
of contamination. Estimates for the possible contamina- 
tion from unresolved point sources are given in Appendix 
\X\ Appendix [Bldetails the creation and processing of the 
gamma-ray sky maps used in this analysis, and provides 
instructions for downloading them. Appendix [Cl contains 
a discussion of Poisson likelihood analysis on smoothed 
maps. 

2. DATA 

The LAT (Large Area Telescope) on Fermi (see 
Geh rels~fc MichelsonI Il999l as well as the Fermi home- 
page^) is a pair-conversion telescope consisting of 16 lay- 
ers of tungsten on top of a calorimeter with a thickness 
of 7 radiation lengths. The entire instrument is wrapped 

* The synchrotron haze depends on the Galactic magnetic field 
while the IC haze depends on the Galactic ISRF, and so the two 
morphologies should not be identical, and could be quite different. 
^ http://ferini.gsfc.nasa.gov/ 



Previous studies of high latitude gamma-ray emis- 
sion have reported measurements of an excess of emis- 
sion abov e the Galactic plane, most notably SAS- 
2 (.Fichtel eTall ligTSl: iKniffen fc Fichtel [l98l. C OS-B 
(iStrond 119841). an d EGRET (iSmialkowski et all Il997t 
iDixon et al.lll998| ). However, in each case the experi- 
ments either covered an insufficient energy range (SAS-2 
and COS-B) or did not have sufficient sensitivity and 
angular resolution (EGRET) to permit a spatial corre- 
lation with the WMAP haze. Fermi overcomes both of 
these obstacles and allows us to search for the gamma-ray 
counterpart to the microwave haze for the first time. 

The presence of an IC signal at the expected level can 
confirm that the microwave haze is indeed synchrotron, 
ruling out the first null hypothesis. From the spectrum 
of the IC, we can estimate the electron spectrum required 
to make the signal, addressing the second null hypoth- 
esis. For example, the presence of ~ 50 GeV IC pho- 
tons requires electrons of > 50 GeV, perhaps much 
greater. Furthermore, IC photons provide valuable in- 
formation about the spatial distribution (disk vs. bulge) 
of the source of these particles, which in turn can con- 
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in a scintillating anti-coincidence detector to provide a 
particle veto. The spacecraft occupies a low Earth orbit 
with an inclination of 25.6°. The field of view is so wide 
that the entire sky may be covered in two orbits by rock- 
ing the spacecraft north of zenith on one orbit and south 
of zenith on the other. Several times per month, Fermi 
interrupts this pattern to point the LAT at a gamma- 
ray burst, though this has little impact on the integrated 
exposure map. When the spacecraft passes through the 
South Atlantic Anomaly (SAA), CR contamination in- 
creases and significant data must be discarded, reducing 
the mean exposure at southern declinations. Beyond a 
zenith angle of 105° the data are significantly contami- 
nated by atmospheric gammas. We excise such data, and 
select only events designated "Class 3" (diffuse class) by 
the LAT pipeline. The LAT collaboration plans to re- 
lease a cleaner class of events in the future, however, at 
the time of this writing, the Class 3 events are the most 
likely to be real diffuse gamma-ray events. 

The events are then binned into a full sky map us- 
ing HEALPix, a convenient iso-latitude equal-area full- 
sky pixelization widely used in the CMB community.^'' 
Spherical harmonic smoothing is straightforward in this 
pixelization, and we smooth the maps to a Gaussian PSF, 
usually of 2° FWHM. The full-sky Fermi maps are dis- 
played in Figure [1] along with an exposure map. See Ap- 
pendix [B] for more details on map construction, smooth- 
ing, masking, and for instructions on how to download 
the maps. 

3. ANALYSIS 

In this work our goal is to test our general preconcep- 
tions about what gamma-ray signals should be present, 
and identify any unexpected features in the Fermi data; 
we avoid detailed comparisons between the data and spe- 
cific theoretical models for the Galactic gamma-ray emis- 
sion. Our approach is to compute linear combinations of 
Fermi maps at several energies and perform template 
analyses with maps of the ISM, radio maps, etc. to see 
what emerges. This sort of open-minded analysis is fiex- 
ible enough to find the unexpected. 

An alternate approach would be to attempt to fit the 
data with a sophisticated physical model, in the context 
of some simulation code (e.g. GALPROP). Such a phys- 
ical model can be quite detailed, including the 3D dis- 
tributions of gas and dust, the 3D distribution of optical 
and FIR photons density and direction, the 3D magnetic 
field, and a 3D model of p, e~ injections. By propagating 
these primary particles with GALPROP, the resulting 
tt", bremsstrahlung and IC signals may be predicted and 
compared with Fermi. However, while a detailed physi- 
cal model will certainly be crucial to a full understanding 
of the Fermi data (and such modeling is currently un- 
derway within the Fermi collaboration), this approach 
may lack the flexibility to identify new emission compo- 
nents that cannot be absorbed by modifying parameters 
in the model. On the other hand, such a model may also 
have too much freedom, so that meaningful patterns are 
absorbed into the flt and left unnoticed. In future work, 
we will take the signals revealed by our initial analysis 

HEALPix software and documentation can be found at 
http://healpix.jpl.nasa.gov, and the IDL routines used in 
this analysis are available as part of the IDLUTILS product at 
http : //sdssSdata. Ibl . gov/sof tware/idlutils. 
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Fig. 2. — GALPROP model illustrating the three primary 
gamma-ray emission mechanisms (see j|3]l and their relative am- 
plitudes in the Galactic plane {\£\ < 30, |fe| < 5). 

and fold them back into a full physical model. 

3.1. Diffuse Gamma Templates 

There are three well-known mechanisms for generating 
gamma-rays at the energies observed by Fermi. First, 
at low (^1 GeV) energies, gamma-ray emission is domi- 
nated by photons produced by the decay of 7r° particles 
generated in the collisions of cosmic ray protons (which 
have been accelerated by SNe) with gas and dust in the 
ISM. Second, relativistic electrons colliding with nuclei 
(mostly protons) in the ISM produce bremsstrahlung 
radiation. Finally, those same electrons interact with 
the interstellar radiation field (ISRF) and inverse Comp- 
ton scatter CMB, infrared, and optical photons up to 
gamma-ray energies. A schematic of the relative impor- 
tance of these emission mechanisms in the Galactic plane 
< 30, 161 < 5) generated by the GALPROP code, 
version 50p (IStrong et al.l l2000t iPorter fc Strond 120051: 
IS trong et al.ll2007( ) is shown in Figured 

Since 7r° gammas and bremsstrahlung are produced 
by interactions of protons and electrons (respectively) 
with the ISM, these emission mechanisms should be mor- 
phologically correlated with other tracers of the ISM, 
such as the SFD du st map based on 100/x far IR data 
(Schl egel et all[l998h . The tt'^ gamma-ray intensity scales 
with the ISM volume density times the proton CR den- 
sity, integrated along the line of sight. In the limit where 
the proton cosmic ray spectrum and density is spatially 
uniform, the ISM column density is a good tracer of 7r° 
emission. Likewise, for a uniform electron spectrum, it is 
a good tracer of bremsstrahlung. Because our analysis is 
limited to |6| > 5°, much of the emission we see is within 
a few kpc, so the assumption of uniform CR density is 
more valid than it would be for the entire Galaxy, partic- 
ularly for protons, which have much larger propagation 
lengths than electrons. 

3.2. Residual Maps 
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Fig. 3. — Residual maps after subtracting the SFD dust template from Fermi maps at various energies. The mask is described in 
Appendix [B] Cross-correlations are done over unmasked pixels and for 75 < £ < 285. Although the template removes much of the 
emission, there is a clear excess towards the Galactic center. This excess also includes a disk-like component which is likely due to IC and 
bremsstrahlung from softer electrons (see Figure [5}. 



Since our goal is to search for an IC emission compo- 
nent with a morphology which roughly matches the mi- 
crowave haze (i.e., centered on the GC, roughly spherical, 
and about 20-40° in radius), we now attempt to remove 
the tt'^ emission from the maps shown in Fig ure [H using 
the sam e ter nplate fitting technique used in iFinkbeined 
(|2004aD and iDobler fc Finkbeineil (|2008[ ). We perform 
multiple types of template fits. Type 1 uses only the 
Fermi map itself at 1-2 GeV which roughly traces 7r° 
emission because the gamma ray sky at those energies is 
dominated by ir^ gammas (with subdominant contribu- 
tions from bremsstrahlung and IC). Type 2 uses only the 
SFD dust map which is a direct tracer of ISM density and 
so should roughly map where the tt*' gammas and much 
of the bremsstrahlung are produced - again, up to some 
uncertainty involving the line of sight distribution. In 
each uniform background is included in the fit, 

making our results insensitive to zero point offsets in the 



maps. 

We model the Fermi map at energy E, F(E), as a lin- 
ear combination of template maps, -Fmodci = Tct, where 



model 



odd 

is a column vector of iVpix unmasked pixel values, 
T is the -/Vtompiato X -^pix template matrix, and the cor- 
relation coefficients ct are chosen to minimize the mean 
squared residual, 

((i^-F„,odoi)')-((F-TcT)2), (1) 

averaging over pixels. The least-squares solution, 

CT{E) = {T^Ty' x{T^F{E)), (2) 

yields the template correlation coefficients at each energy. 
In this fit, we mask out the Fermi 3- month point source 
catalog as well as the LMC, SMC, Orion-Barnard's Loop, 
and Cen A (see Appendix iBjfor details) . We also mask all 
pixels with Galactic latitude |6| < 5°. Cross-correlations 
are done over unmasked pixels and for several different 
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Fig. 4. — Top: the WMAP Q-band (41 GHz) total synchrotron (left) and haze (right). Bottom: the residual map obtained with the 
Type 1 template fit using the SFD map (left) and the 1-2 GeV Fermi map (right) as a template for tt^ decay emission. These residual 
maps are the same as the right colum n second r ow panel of Figure [3] and Figure |5] respectively though with a different stretch. Both maps 
are shown with the mask used in the Dobler & : Finkbcincr ( 2008) microwave analysis for comparison. The upper panel represents cosmic 
ray electrons interacting with the Galactic magnetic field to produce synchrotron while the lower panel represents cosmic ray electrons 
interacting with the ISRF to produce IC emission, and colliding with the ISM to produce bremsstrahlung. The morphological similarities 
between the microwave haze and gamma-ray haze (right) are striking, indicating a correspondence between the two electron populations. 



ranges in i: eight longitudinal slices that have \£\ > 75° 
(to avoid the GC) and width 30° (Regions 1-7) as weU 
as one region with 75° < £ < 285° (Region 8). Note that 
Region 8 is the union of Regions 1-7; Regions 1-7 are fit 
individually to show the variation in the fit spectrum, 
and Region 8 is fit to obtain the mean outer Galaxy 
signal. With these correlation coefficients, we define the 
residual map to be, 

RriE) = F{E) - ct{E) x T. (3) 
To the extent that the templates in T match the morphol- 
ogy of the 7r° and bremsstrahlung gammas, the residual 
map will include only IC emission. 

Figure [3] shows the resultant residual maps using the 
SFD map as a morphological tracer of 7r° emission for 
the Region 8 fits. The most striking feature of the dif- 
ference maps is the extended emission centered around 
the Galactic center and extending roughly 40° in b. 
The morphological correlation between the WMAP syn- 
chrotron and the i?sFD(5 — 10 GeV) is striking as is 
shown in Figure S) Here the 41 GH z synchrotron (haze 
plus Haslam-correlated eni ission, see lHaslam et al.|[l982l : 
iDobler fc Finkbeineill2008[ ) is shown side by side with the 
5-10 GeV F ermi residual map with the mask used in the 
iDobler fc~F inkbcincr (2008) microwave analysis overlaid 
for visualization. 

Figure [5] shows residual maps using the 1-2 GeV Fermi 
maps as a template for the tt" emission. This sort of 
"internal linear combination" has the advantage that 7r° 
emission cancels out as long as the shape of the proton 
CR spectrum is the same everywhere - it does not rely on 
the proton CR density to be uniform. The residual maps 
look largely similar to the case with the SFD template 
regressed out, but there are some notable differences. In 
particular, using this 1-2 GeV template, the IC haze has 



a slightly "taller" appearance. This seems to be due to a 
disk-like component that is present in the 1-2 GeV maps 
but not in the dust map. This is probably because the 
1-2 GeV Fermi map is not entirely tt" emission, but also 
contains bremsstrahlung and IC components generated 
by electrons accelerated by SNe in the disk. The re- 
sult is that when this lower energy map is regressed out 
(i.e., cross-correlated and subtracted) from higher energy 
maps, this emission component is subtracted along with 
the dust-correlated emission. Conversely, the SFD map 
contains only emission from dust grains and not rela- 
tivistic electrons; while the bremsstrahlung from those 
electrons largely traces the gas distribution, the IC does 
not, and so is not regressed out when using the SFD 
template. 

Despite the significant shot noise in the Fermi map. 
Figure |4] shows that there is a clear morphological corre- 
lation between the microwave haze and the gamma-ray 
haze. Of course, we do not expect the morphologies to 
agree perfectly since the microwave haze is generated by 
interactions of electrons with the magnetic field while this 
IC haze is due to interactions of the electrons with the 
ISRF. Nevertheless, this is evidence that the microwave 
haze seen in the WMAP data is indeed synchrotron and 
is not some other component such as free-free emission 
or spinning dust. 

3.3. Four- Component Template Fits 

Because the solution to Equation [5] minimizes the vari- 
ance of the (zero mean) residual map R{E), the pres- 
ence of the IC haze affects the coefficients ct{E), since 
there is non-zero spatial correlation between the tem- 
plates and the IC haze. To relax the stress on the fit, we 
expand our template analysis with a third Type. Type 
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Fig. 5. — The same as Figure[3]but using the Fermi 1-2 GeV map for cross-correlations instead. Unlike the SFD dust map which should 
trace tt" (and subdominant bremsstrahlung) emission only, the low energy Fermi map includes the soft IC associated with lower energy 
electrons. In fact comparing the residuals in this figure with those in Figure (3] it is clear that the disk-like component has been subtracted 
leaving only the IC haze. Furthermore, the IC haze is more prominent in the high energy maps indicating a harder spectrum than tt'' 
emission which is the dominant emission mechanism at ~1 GeV energies. 



3 uses four templates in T: the SFD ma p to trace 7r° 
and bremsstrahlung emission, the 408 MHz'Hasla m et al\ 
([1982) map which is dominated by radio wavelength syn- 
chrotron and thus roughly traces soft spectrum electrons 
which produce soft IC and bremsstrahlung, and a bi- 
variate Gaussian of width ae — 15° and ai, — 25°. We 
note that this template is chosen to roughly match the 
morphology in Figure [5] and has no other physical mo- 
tivation. We also use a uniform template to fit out the 
isotropic background signal in the maps, again, making 
our results insensitive to zero points. Lastly, for this fit 
we use all values in £ (Region 9). 

Note that since the bremsstrahlung originates from in- 
teractions of the electrons with the ISM, its spatial dis- 
tribution depends on both the gas density and the cos- 
mic ray electron density; consequently, some contribution 
from bremsstrahlung will be present in both the SFD- 
correlated and Haslam-correlated emission. 



The previous fits were done with uniform weighting 
and assuming Gaussian errors, minimizing x^- For the 
Type 3 fit we do a more careful regression, maximizing 
the Poisson likelihood of the 4-template model in order to 
weight the Fermi data properly. In other words, for each 
set of model parameters, we compute the log likelihood 



ln£ = ^ fcj In/ij - - \n{kil), 



(4) 



where /ij is the synthetic map (i.e., linear combination 
of templates) at pixel i, and k is the map of observed 
counts. Note that the last term does not depend on the 
model parameters. It may appear strange at first to com- 
pute a Poisson likelihood on smoothed maps, however, 
the smoothing is necessary to match PSFs at different 
energies and with various templates (some of which have 
lower resolution than Fermi in the energy range of inter- 
est). The smoothing itself does not pose any problems 
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Fig. 6. — The templates and fit solutions used in the Type 3 (see i]3.2ll template fits. Upper left: the SFD dust map, upper right: 
the Haslam 408 MHz map, middle left: the bivariate Gaussian haze template, middle right: the Fermi map at 10-20 GeV (same as the 
first column, third row of Figure [5] but with a different stretch to show the detailed morphological structure), lower left: the best fit 
template solution for the observed emission, lower right: the residual map. Note the very small residuals indicating that the template fit 
is a remarkably good representation of the data over large areas of the sky. 



for relative likelihoods, as we show in Appendix [Cl 

However, we must keep in mind that the uncertainties 
derived in this way are the formal errors corresponding 
to Aln£ ~ 1/2, which would be 1 cr in the case of Gaus- 
sian errors. The error bars plotted are simply the square 
root of the diagonals of the covariance matrix. This es- 
timate of the uncertainty should be accurate at high en- 
ergies, where photon Poisson noise dominates. At low 
energies, although the formal errors properly reflect the 
uncertainty in the fit coefficients for this simple model, 
the true uncertainty is dominated by the fact that the 
4-template model is not an adequate representation of 
the data. 

Figure |6] shows the skymaps and best fit solution in- 
cluding the residual map at 10-20 GeV while Figure [7] 
shows residual maps at other energies. It is clear from 
these residuals that the template fitting produces a rel- 
atively good approximation of the gamma-ray data over 
large areas of sky. Furthermore, Figure [7] shows that 
not including the bivariate Gaussian template for the IC 
haze yields a statistically significant residual towards the 
center indicating that a model including an IC haze is a 
better match to the data then one without. The promi- 
nent North Polar spur feature in the Haslam map, which 



is thought to origi nate from synchro tron emission from 
electrons in Loop I (jLarge et al.lll962l l , is over-subtracted 
in each case, because the North Polar spur is brighter in 
the Haslam map than in the gamma-ray maps (i.e. the 
ratio of synchrotron microwaves to IC gamma rays in the 
North Polar spur is larger than in the rest of the Haslam 
map)." This may be due to different ISRF, B-field and 
ISM density values in Loop I relative to the inner Galaxy 
(since Loop I is thought to be quite nearby), or may be 
due to a softer-than-usual electron spectrum in Loop I, 
since the electrons producing the synchrotron measured 
in the Haslam map are much lower energy than those pro- 
ducing IC gamma-rays at energies measured by Fermi. 

Figure [8] shows these same residual maps, but with a 
smoothing of 10° which is on the order of the scale of the 
haze emission. With this large smoothing, smaller scale 
variations (due to individual photons at high energies) 
are smoothed over and the residual maps clearly show 
that the haze is a robust feature at all energies. 

In Figure [9] we show the results of a four template fit 
using the 1-2 GeV map instead of SFD to trace the tt" 



See ICasandiian et al.l II2009I) for a discussion of gamma-rays 
from Loop I seen by Fermi. 
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emission. This figure shows that the haze is not due to 
the SFD template being an imperfect tracer of emis- 
sion. If the proton cosmic ray density is higher towards 
the GC, then SFD may systematicahy underestimate the 
tt" emission there and perhaps the gamma-ray haze is the 
resuh. However, Figure |9] shows that this is clearly not 
the case. The 1-2 GeV template includes the effects of 
proton cosmic ray density variations (as well as line of 
sight gas density effects) and the haze remains as a ro- 
bust residual. That is, the haze is not due to imper fect 
templates as sugg ested bv lLinden fc Profumol (|2010[ ). 

In Figure [10] we show the haze amplitude (residual 
maps from the Type 3 fit plus the correlation-coefRcient- 
weighted bivariate Gaussian template) as a function of 
Galactic latitude and for the longitudinal bin: —15° < 
^ < 15°. Although the data are noisy, the figure shows 
that the Fermi haze dies off by roughly h ^ 50° in all 
energy bands. For comparison, we show the same plot 
for just the residual map. The figure also shows the am- 
plitude as a function of longitude for the latitudinal bin 
-20° <h < -10°. The more rapid fall off of the haze 
emission with I compared to h indicates a haze morphol- 
ogy elongated in the b direction. 

There are two important features to note about Figures 
[3] through [71 First, the IC haze has a spectrum which is 
harder than the other IC in the Galaxy. This is evidenced 
by the fact that using the 1-2 GeV Fermi map as a tem- 
plate removes the IC emission from disk electrons but the 
haze IC fades more slowly with energy. Second, the IC 
haze is morphologically distinct from either the tt*^ emis- 
sion or IC and bremsstrahlung from the disk electrons. 
These two facts taken together strongly suggested that 
the electrons responsible for the microwave and gamma- 
ray haze are from a separate component with a harder 
spectrum than SN shock-accelerated electrons. 

3.4. Fermi Galactic Diffuse Emission Model 

So far, we have only done template fits of maps of 
data to other maps of data. Using the SFD dust map as 
a tracer of 7r° emission is equivalent to assuming that the 
proton CR spectrum and density are spatially uniform, 
and the dust/gas ratio is constant. Using the Haslam 
map to estimate IC emission is equivalent to assuming 
that the B-field and ISRF have similar spatial varia- 
tion, and neglecting the anisotropies in both IC and syn- 
chrotr on emission. Consequen tly, there have been con- 
cerns ([Linden fc Profumoll2010[ ) that this template-based 
analysis might introduce spurious large-scale residuals. 
To address these concerns, we investigate whether the 
Fermi diffuse model contains a structure similar to the 
haze emission. 

The Fermi team provides a model of diffuse gamma- 
ray emission consisting of maps sampled at 30 energy 
bins from 50 MeV to 100 GeV.^^ These maps are based 
on template fits to the gamma-ray data and also include 
an IC component generated by the GALPROP cosmic 
ray propagation code. 

The difference between the model and data is shown 
in Figure [TlJ We have interpolated the model to the 
energy ranges of interest and performed the simple one- 
template fit to the data (analogous to our Type 1 and 2 

The background models can be downloaded from 
http : //f ermi . gsf c .nasa. gov/ssc/data/access/lat. 



fits described above). This fit allows for any error in the 
normalization of the diffuse model; the fit coefficients are 
within a few percent of unity in every case. 

At each energy range, the haze is clearly visible in the 
residual maps. GALPROP uses the standard inhomo- 
geneous ISRF model (Por ter fc Strong 2005). making it 
unlikely that the observed residual is due to the expected 
spatial variation of the ISRF. Furthermore, the diffuse 
model includes a disk-like injection of primary electrons 
and estimates for the line of sight density variations of 
the interstellar medium. The fact that the haze residual 
remains suggests that, given the propagation model in- 
cluded in GALPROP, disk-like sources will not produce 
the observed IC emission, nor will expected variations in 
the gas density and proton CR density along the line of 
sight. 

3.5. Spectra 

While the morphology of the gamma-ray haze is indica- 
tive of IC emission from the microwave haze electrons, we 
now attempt to estimate the spectrum of this emission. 
This is difficult for several reasons. First, tt*^ emission 
is dominant (or nearly so) at most energies in FermVs 
energy range. Thus, in a given region, we must estimate 
the spectral shape and amplitude of the i:^ emission in 
order to subtract it from the total. Second, the total 
number of photons measured by Fermi decreases rapidly 
with increasing energy. For example, in the inner Galaxy 
< 30, |6| < 5), there are only ^4300 photons between 
10 and 100 GeV, resulting in substantial Poisson errors. 
Lastly, there is also significant isotropic background due 
to extra-galactic emission and charged particle contami- 
nation, including heavy nuclei at high energies. We will 
show that the isotropic backgrounds are manageable, and 
present two types of spectra: template coefficients ct{E) 
and total fiuxes dN/dE. 

3.6. Isotropic Background 

The Fermi data contain gamma-rays from an un- 
resolved extragal actic signal with dN/dE ~ i?"^-^^ 
(|Ackermannll2009| ) as well as particle contamination. We 
make no attempt to separate the extragalactic gamma- 
ray signal from the particle contamination. Instead we 
measure the spectrum of this (nearly) isotropic back- 
ground in 8 regions at high latitude and test the assump- 
tion of isotropy. Specifically, we take combinations of 
longitude ranges of -180° < £ < -90°, -90° < ^ < 0°, 
0° <e < 90° and 90° < ^ < 180° together with latitude 
ranges of -90° < 6 < -60° and 60° < 6 < 90° for the 
4x2 = 8 regions. The regions are chosen to be far from 
the plane to avoid contamination from Galactic emis- 
sions. We use the point-source masked maps described 
above, binned in 12 energy bins from 0.3 - 300 GeV. 
The energy bins, mean background, etc. are given in Ta- 
ble [T] and the 8 spectra, along with the inverse- variance 
weighted mean, are plotted in Figure [12j We note that 
the error bars are the standard deviation of the 8 regions, 
not the standard deviation of the mean. To be conserva- 
tive, we use this standard deviation as the uncertainty in 
the background in the remaining stages of the analysis. 
Future work by the Fermi team to understand both the 
particle contamination and the gamma-ray sky will likely 
reduce this error in the background substantially. 
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Fig. 7. — Residual maps using the Type 3 template fit. The right column is the same as the lower right map in Figure [6] but for maps 
at different energy bands. The left column performs the same fit without including a bivariate Gaussian template for the IC haze. It is 
clear that not including the haze template results in a significant residual towards the GC in each energy band, but particularly at high 
energies. Including the haze template improves ln£ by 504, 215, 78, and 54, respectively, for the 4 energy bins shown. 



3.7. Template- Correlated Spectra 

Figure [T^ shows ct{E) x (T) for the two templates and 
regions 1-7 used in the Type 1 and Type 2 fits along with 
the model 7r° s pectrum from GALPROP, which uses the 
iBlattnig et al.l (2000) paranicterizations for pion produc- 
tion. It is clear from the figure that the cross-correlation 
technique produces tt^ spectra that are remarkably simi- 
lar to the model spectrum at low energies, while at high 
energies the cross-correlation spectrum is slightly higher 
than the model spectrum. This could be due to a number 
of reasons such as non-zero spatial correlation between 
the templates and the harder spectrum haze IC, contam- 
ination from background events like heavy nuclei, or un- 
certainties in the 7r° emission model. Of these, the first is 
most likely since the cross-correlation between the tem- 
plates and a nearly isotropic background is likely small 
and since the spectrum of 7r° gammas is quite well known. 

Template-correlated spectra for the Type 3 template 



fit are shown in Figure [M) Here the correlation coeffi- 
cients are weighted by the mean of each template in the 
"haze" region (see Table [2]). As shown in the figure, the 
spectra for the SFD and Haslam maps reasonably match 
the model expectations in that region^'^. That is, the 
SFD-correlated emission roughly follows the model 7r° 
spectrum while the Haslam-correlated spectrum resem- 
bles a combination of IC and bremsstrahlung emission. 
However, the /laze-correlated emission is clearly signif- 
icantly harder than either of these components. This 
fact coupled with the distinct spatial morphology of the 
haze indicates that the IC haze is generated by a separate 
electron component. 

3.8. Total Intensity Spectra 

The GALPROP model here was tuned to match locally mea- 
sured protons and anti-protons as well as locally measured electrons 
at ~ 20-30 GeV. 
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Fig. 8. — The same as Figure[7]but for a FWHM smoothing of 10°. Smoothing the residuals at this scale demonstrates that the haze 
is a robust feature and, in particular, is not the consequence of single photon fluctuations at high energies. 



While the template-correlated spectra and residual 
maps are useful for identifying separate components, for 
the purposes of comparing the map intensities to a model 
for the physical mechanisms, we now generate total in- 
tensity spectra in several regions of interest. We define 
three key regions: a "haze re gion" south of the GC^'', 
the "four corners" region from lCholis et all (I2b 09b'). and 
a Galactic plane region used bv iPorte r (2009). 

The spectrum of each region is shown (Figure I15p . 
with the background spectrum from Figure [T2l subtracted 
from the other spectra to remove any isotropic compo- 
nent. The inner Galaxy region clearly shows the low en- 

The region south of the GC is greatest interest for studying the 
microwave and gamma-ray haze. This is because, north of the GC, 
the microwave maps from WMAP include bright free-free emission 
from f Oph, and spinning dust emission from p Oph, both of which 
are bright enough to leave substantial residuals, even though they 
are relatively well subtracted by the template fitting described in 
IDob lcr & Finkbciner (2008). South of the Galactic center there are 
some small dust and gas features that should provide some signal 
in the gamma-ray map, but the situation is much simpler. 



ergy behavior characteristic of tt" emission (though it is 
important to note that there may be significant emission 
from unresolved point sources such as pulsars in this re- 
gion wh ich have a similar spectral shape; see Appendix 
Eland lAbdo et al.ll2009cD . At high energies however, 
there is a significant excess which is not expected for tt" 
emission (cf., Figure [2]). The spectrum of this excess is 
similar to the spectrum derived for the haze template 
from Figure [TJ] and is consistent with an IC signal from 
a hard electron population with energies > 10 GeV. The 
evidence for this excess is more pronounced in the haze 
and four corners region. Lastly, we take the template es- 
timate of TT^ emission in the haze region from the Type 3 
fits and subtract it from the total emission. This reveals 
two clear features: a bump centered on roughly 1-2 GeV 
that is likely due to either bremsstrahlung from a low 
energy electron component or emission from unresolved 
pulsars, and a hard tail above ^ 10 GeV that is the IC 
signal from the haze electrons. 
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Fig. 9. — The same as Figure[7]but using the 1-2 GeV map instead of the SFD dust map to trace tt" emission. The clear haze residual 
seen, particularly at high energies, indicates that the gamma-ray haze is not due to shortcomings of the SFD template resulting from 
variations in the proton cosmic ray density. These variations, as well as line of sight gas density effects, are automatically included in the 
1-2 GeV template. 



3.9. Comparison to the Microwave Haze 

Both the morphology and the relatively hard spectruni 
of the gamma-ray haze motivate a common physical ori- 
gin with the WMAP haze. We now provide a simple 
estimate of the microwave and gamma-ray signals from 
a population of hard electrons in the inner Galaxy, to 
demonstrate that the magnitudes and spectral indices of 
the two signals are consistent for reasonable parameter 
values. 

We consider a steady-state electron spectrum described 
by a power law, dN/dE oc E^", with a high-energy cutoff 
at 1 TeV (here the cutoff is implemented as a step func- 
tion, not an exponential fall-off; of course this is only 
an approximation to the true spectrum). This choice 
is motivated by the local m easurement of the c osmic ray 
electron spectrum by Fermi (jAbdo et al.ll2009"bl ) . We con- 
sider a region ^ 2 kpc above the Galactic center, as an 
example (and since both hazes are reasonably well mea- 



sured there), and employ the model for the ISRF used in 
GALPROP (jPorter fc Strondl2005f ) at 2 kpc above the 
center. We normalize the synchrotron to the approxi- 
mate value measure d by WMAP in the 23 GHz K-band 
(jHooper et al.l [20071 ) . 15° below the Galactic plane, 
and compute the corresponding synchrotron and IC spec- 
tra. The WMAP Haze was estima ted to have a spec- 
trum U oc 15 = 0.39 - 0.67 ( Dobler k Finkbeined 
l2008f l. corresponding approximately to an electron spec- 
tral index of a ~ 1.8 — 2.4; Figure [16] shows our re- 
sults for a magnetic field of 10 /iG, and electron spectral 
indices a = 2 — 3. This field strength is appropriate 
for an exponential model for the Galactic S-field inten- 
sity, \B\ = \Bo\e-^/^% with 5n » 20 - 30/^G (which 
seems reasonable; see e.g. iFerrierd (|2009D and references 
therein) and scale height Zg ~ 2 kpc (the default value 
in GALPROP). We find good agreement in the case of 
a w 2 — 2.5, consistent with the spectrum of the WMAP 
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Fig. 10. — Top left: The haze amplitude (Type 3 residual map plus haze template) as a function of Galactic latitude for four different 
energy ranges. The data are binned in steps of 0.055 in sin 6, i.e. roughly 4° bins near the plane. The plot shows that haze amplitude has 
roughly the same fall off with b (~ 40°) in all energy ranges. Top right: The same but for the Type 3 residual map. Bottom left and right: 
the same as the top left and right, but as a function of I for a fixed range of b. The haze profile clearly falls off more quickly with I than 
with b indicating a profile elongated in the b direction. 



Haze. 

3.10. Comments on Haze Morphology 

Although a detailed analysis of the possible sources of 
the Fermi haze are beyond the scope of this paper, a few 
simple comments are in order. First, the profile is not 
well described by a disk source. While quantifying this 
is a subtle task, the success of the template makes this 
point clear. There are many possibilities to explain the 
oblong shape, should that persist in future data. For in- 
stance, AGN jets and triaxial DM profiles could both pro- 
duce signals of this shape. Even approximately spherical 
injection could yield such a signal, should the diffusion 
be considerably anisotropic. That said, while one might 
have attempted to invoke, e.g., rapid and significant lon- 
gitudinal variation of the magnetic field to explain the 



microwave haze with a disk-like electron injection pro- 
file, such approaches are no longer tenable (and moreover 
already had significant tension with an understanding of 
the Haslam synchrotron maps). The presence of this fea- 
ture in both gamma rays and microwaves demonstrates 
that the electrons themselves do not follow a disk-like 
injection profile. Furthermore, we emphasize that the 
angular scale of the haze is very large, roughly 25-40° in 
h. This region of the sky is far from the Galactic disk or 
the very GC where complications from the central black 
hole, other point sources, or disk-like pulsars could com- 
plicate the analysis. 

4. DISCUSSION/CONCLUSIONS 

We have presented full sky maps generated from pho- 
ton events in the first year data release of the Fermi 
Gamma-Ray Space Telescope (see Appendix |B] for data 
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Fig. 11. — The same as Figure [3] but using the official Fermi team Galactic diffuse model. The strong haze residual in the right hand 
panels shows that the haze is not included in the model. Since the IC emission in the model was obtained with GALPROP, this indicates 
that variation in the ISRF with position in the Galaxy cannot account for the haze emission. In addition, the haze morphology is not 
reproduced by disk-like injection of electrons, nor by the cosmic ray propagation model employed by GALPROP. 
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TABLE 1 



E range 


Energy 


background 




(GeV) 


(GeV) 


(xlO °GeV/cm^/s/sr) 




0.3 — 0.5 


0.4 


2 235 + 1 87 


86.467 


0.5 - 0.9 


0.7 


1.901 ± 0.190 


76.812 


0.9 - 1.7 


1.3 


1.657 ±0.181 


51.885 


1.7 - 3.0 


2.2 


1.340 ±0.157 


27.169 


3.0 - 5.3 


4.0 


1.093 ± 0.119 


10.981 


5.3 - 9.5 


7.1 


0.886 ± 0.069 


2.720 


9.5- 16.9 


12.7 


0.687 ±0.098 


4.190 


16.9 - 30.0 


22.5 


0.632 ± 0.068 


1.357 


30.0 - 53.3 


40.0 


0.695 ± 0.078 


0.936 


53.3 - 94.9 


71.1 


0.624 ±0.125 


1.183 


94.9 - 168.7 


126.5 


0.704 ±0.199 


1.594 


168.7 - 300.0 


225.0 


0.420 ±0.205 


1.949 



Note. — The background is tabulated in each energy bin 
for 8 polar (|b| > 60°) regions after point source masking. The 
uncertainty is the RMS of the 8 regions, not the uncertainty 
in the mean, which is \/7 times smaller. For each region 
we also compute the Poisson uncertainty, o-p, which for low 
E is much smaller. To test the significance of background 
variation, we compute xt = {'^'^ 1'^'^)^ averaging over the 8 
regions. There is no strong indication of anisotropy at high 
latitude for £ > 5 GeV. 
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Fig. 12. — Background spectrum, including gamma-rays and 
particle backgrounds. The background is tabulated in each energy 
bin for 8 high-latitude polar (|fe| > 60°) regions after point source 
removal (dotted gray lines). The mean of the 8 regions is shown 
with a solid black line and the uncertainty is the RMS of the 8 re- 
gions, not the uncertainty in the mean, which is \/7 times smaller. 
For reference, a dN/dE oc E~'^-^^ model for the unresolved back- 
ground is shown with a dot-dashed line. 



TABLE 2 



Region 


£ range 


b range 


Inner Galaxy 
Haze 
4 Corners 


\e\ < 30° 
\e\ < 15° 
5 < |£| < 10° 


\b\ < 5° 
-30 <b< -10° 
5 < |6| < 10° 



Note. — Regions for which the total gamma- 
ray intensity is evaluated (see also Figure [TSl l. 



processing details). 

Using a template fitting technique, we have approx- 
imated both the spectrum and morphology of the well 
known gamma-ray emission components at Fermi en- 
ergies. The SFD dust map was used to trace the tt*^ 
decay gammas generated by collisions of cosmic ray pro- 
tons with the ISM, while the Haslam 408 MHz map was 
used to trace inverse Compton (IC) scattered photons 
from interactions of supernova shock-accelerated elec- 
trons (~ 1 — 10 GeV) with the interstellar radiation field 
(ISRF). Bremsstrahlung radiation, generated by interac- 
tions of these electrons with the ISM, should be approx- 
imately traced by some combination of these two maps. 
Although our template fitting technique is subject to sig- 
nificant uncertainties due to uncertain line of sight gas 
and CR distributions, a robust positive residual has been 
identified. 

This excess diffuse emission is centered on the Galac- 
tic center, and can be parameterized by a simple two- 
dimensional Gaussian template {at = 15°, ab = 25°). 
The template-correlated spectrum of this emission is sig- 
nificantly harder than either tt'^ emission or IC from 
softer electrons, whose fitted spectra agree well with 
models. This harder spectrum coupled with the distinct 
spatial morphology of the gamma-ray and microwave 
haze are evidence that these electrons originate from a 
separate component than the softer SN shock-accelerated 
electrons. 

The gamma-ray excess is almost certainly the IC 
counterpart of the microwave haze excess described by 
[Unkbcineil (|2004aD and Doblcr & Finkbeincr (2008). Al- 
though it is still possible that a significant fraction is 
prompt photons from W IMP annihilations (e.g. the 200 
GeV wino advocated bv lGraiek et al.ll2009f ) such expla- 
nations are difficult to reconcile with the spatial simi- 
larity to the WMAP haze (see Figure |4]). The simplest 
hypothesis is that the signal is mainly IC from the same 
electrons that produce the WMAP haze synchrotron. 

This addresses the stubborn question about the ori- 
gin of the WMAP haze. Until recently, it has been ar- 
gued that the WMAP haze had alternative explanations, 
such as free-free emission from hot gas or spinning dipole 
emission from rapidly rotating dust grains. However, the 
existence of this IC signal proves that the microwave haze 
is indeed synchrotron emission from a hard electron spec- 
trum. 

Fermi LAT photon data are contaminated by particle 
events, especially at high energies. We have taken care to 
account for the isotropic background resulting from ex- 
tragalactic sources, cosmic ray contamination, and heavy 
nuclei contamination and found that this background, 
though significant, is below the observed IC excess even 
up to 100 GeV. Particle contamination is extremely un- 
likely to mimic the observed signal. 

The LAT collaboration continues to refine the cuts 
used to define "diffuse class" events, and plans to re- 
lease a cleaner class of events in coming months. This, 
along with a new public version of GALPROP, including 
updated ISRF models, will allow a more sophisticated 
analysis than that presented in this paper. We eagerly 
await the release of these software and data products. 

The spectrum and morphology of both the microwave 
and gamma-ray haze constrain explanations for the 
source of these electrons. There have been speculations 




Fig. 13. — Left: The tt" spectrum derived using the cross-correlation technique defined in i|3.2l and with the Fermi 1-2 GeV map as a 
tt" template. The dotted Gray lines are the cross-correlation spectra in several diff'erent regions of the sky and the solid black line is the 
moan of those spectra. Error bars on the cross-correlation spectrum are defined as the variance in the values for the different regions. The 
dot-dashed line is the tt'^ spectrum output from a GALPROP model (shown here with arbitrary normalization). Right: The same but 
using the SFD dust map as a template. 

that the microwave haze could indicate new physics, such 
as the decay or annihilation of dark matter, or new as- 
trophysics, such as a GRB explosion, an AGN jet, or a 
spheroidal population of pulsars emitting hard electrons. 
We do not speculate in this paper on the origin of the 
haze electrons, other than to make the general observa- 
tion that the roughly spherical morphology of the haze 
makes it difficult to explain with any population of disk 
objects, such as pulsars. The search for new physics - or 
an improved understanding of conventional astrophysics 
- will be the topic of future work. 
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Fig. 14. — Correlation coefficients for the templates used in the Type 3 template fit (see il3.2p . Upper left: SFD-correlated spectrum 
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Lower left: haze template-correlated emission. This component has a notably harder spectrum than both the SFD- and Haslam-correlated 
spectra or their model shapes (dashed lines, cf., Figure indicating a separate component. The dashed lines in the upper two panels are 
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APPENDIX 
UNRESOLVED POINT SOURCES 

In this section we explore the possibihty that the diffuse excess discussed in this work (the Fermi haze) could 
originate from a large number of unresolved point sources. The limit where the contribution from unresolved point 
sources is dominated by emission from many very faint 1 count) sources cannot be distinguished from purely diffuse 
emission; however, if a smaller number of brighter 1 count) point sources dominate, photon events coincide within 
the point spread function (PSF) at a greater rate than would be the case for purely diffuse emission. 

We model the luminosity function of unresolved point sources by a simple power law for the expected number of 
counts, S: dN/dS cx with S'min < -S* < Smax- Smaller values of a, and larger values of S'max and 5111111, correspond 
to a larger fraction of bright point sources and thus strengthen any upper bounds on the flux fraction from point 
sources. Examining known point source populations (of AGN and X-ray binaries) generally yields spectral indices 
a '--^ 1.5 — 2.2, although at high luminosities steeper power laws with a ^ 2.5 have been observed in flat spectrum radio 
quasars, and at low luminosities very shallow power laws with a ^ 1.2 have been measured in LMXB populations 
(IGrimm et al.l[200a iGilfMlovl [20Q1 IKim fc Fabbicmol [200l IVoss fc Gilfanov! [20061: IPadovani et all [20071: lAbdo etHI 
l2009al) . We take a = 2.2 and S'max = 10, as pessimistic benchmark parameters to provide a robust upper bound. 

This power law must break at some low- luminosity cutoff to avoid divergences if a > 2; we will display the statistical 
limits obtained for Smin = 0.1 and 1 count. It is possible that any remaining excess emission could originate from 
point sources with expected counts/year below 0.1, but this would require the number of point sources to exceed the 
number of unexplained excess photons by greater than a factor of 10. 

To bound the fractional flux originating from point sources, we apply a statistical test of isotropy described in detail 
bv lSlatver fc Finkbeineii (|2009[ ). designed for the case where the density of events is < 1 count per PSF, to the Fermi 
LAT data in the haze region < 15, —30 < b < —10). In brief, we compute (1) the fraction of "isolated events" 
with no neighbors within some test radius r, and (2) the fraction of 10^ randomly distributed circles of radius r which 
contain no events ("empty circles"). If r is chosen appropriately (i.e. r ^ PSF) then the ratio of these two quantities, 
R = "isolated" / "empty," is related to the fraction of counts arising from unresolved point sources. If the average 
number of events per circle of radius r, denoted A, is greater than 1, then we redefine "isolated" events to be those 
with fewer than A neighbors, and "empty circles" to be those containing fewer than A events. Neighbors are found 
efficiently using the publicly available IDL routine spherematch, and we take the test radius r to be equal to the 
estimated la value for the LAT PSF, corresponding to 39% flux containment^^ (our results are not very sensitive to 
this choice). We apply this test separately to the diffuse class data in each of the energy bins in Fig. [1] removing events 
with a large zenith angle as described previously, and calibrate our results via Monte Carlo simulations. The smaller 
the PSF, the better the bounds (the PSF is an input to our Monte Carlo calibration); since the PSF size decreases 
with increasing energy, and is greater for back-converting events compared to front-converting events, for each energy 
bin we take the PSF for back-converting events at the minimum energy of the bin, in order to set robust limits. 



The energy-dependent PSF for the LAT is taken from http : //www-glast . slac . Stanford. edu/software/IS/glast_lat_performaiice .htm. 
The radius of 68% containment used by the Fermi Collaboration would in the case of a Gaussian PSF be rgg = 1.51a. 
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Fig. 17. — Upper limits on the fraction of diffuse class events originating from unresolved point sources, in the Fermi haze region 
(|Z| < 15, —30 < 6 < —10), at 95% confidence. The dark shaded region indicates bounds on sources with an expected contribution of 1+ 
counts/year, while the light shaded region includes sources with 0.1+ counts/year; both cases assume a luminosity function with spectral 
index a = 2.2 and 5max = 10. The dotted and dashed lines indicate the bounds with a steeper luminosity function a = 1.8, for Smin = 1 
and 5niin = 0.1 res pect ively. Data points indicate the fraction of emission correlated with the "Fermi Haze", obtained from the spectra 
summarized in Fig. 1141 The table shows the value of Icr for the detector point spread function used in the Monte Carlo calibration, as 
conservatively estimated for each energy bin. 

At low energies we find values of R consistent with a significant point source component (although the point source 
fraction cannot be precisely determined without knowing the luminosity function). Above 10 GeV, however, the value 
of R is consistent with entirely diffuse emission. In the 50 — 100 GeV bin, low statistics nonetheless allow the haze- 
correlated fraction of the emission to be explained entirely by point sources with 0.1+ expected counts/year (since in 
this bin the haze-correlated flux in this region corresponds to ^ 30 photons), but at ^ 10 — 50 GeV there is a clear 
excess over the 95% confidence limits on the unresolved point source contribution, for S'min = 0.1, 1 (Fig. [T7)) . There 
are 900 diffuse class events in this energy range and spatial region which pass the zenith angle cut, of which ^1/3 are 
attributed to the haze: in order for point sources providing an expected flux of less than 0.1 counts/year to make up 
the difference, there would have to be 1000-1- such sources emitting 10-50 GeV gamma rays in this 600-square-degree 
region of the sky (in addition to similar source populations in the other regions where the haze-correlated emission is 
bright). Therefore, even for these conservative assumptions, it seems very unlikely that the hard spectral shape of the 
Fermi haze could be due to point source contamination over the entire relevant energy range. 

FERMI FULL-SKY MAP PROCESSING AND DATA RELEASE 

In this section we describe full-sky maps generated from photon events in the first year data release of the Fermi 
Gamma-Ray Space Telescope. These maps are corrected for exposure, point source masked, and are smoothed to a 
Gaussian PSF, usually of 2° FWHM. The lowest energy maps are smoothed to 3° or 4°. For convenience, we provide 
two sets of energy bins: 12 logarithmically spaced bins for use in spectral analyses ("specbin") and 8 somewhat larger 
bins with better signal/noise for visual inspection ("imbin"). 

Although the maps were generated from a public Fermi data release, the maps are not an official LAT data release, 
and the procedure used to make them has not been endorsed by the LAT collaboration. 



Event Selection and Binning 

The LAT is a pair-conversion telescope, in which incoming photons strike layers of tungsten and convert to e+e~ 
pairs, which are then tracked (to determine direction) on the way to a calorimeter (to determine energy). The first 
year P6_V3_DIFFUSE data file contains records of 15,878,650 events, providing the time, energy, arrival direction, (with 
respect to the spacecraft, and in celestial coordinates), and zenith angle for each event. Events are divided into 3 
classes, using a number of cuts to separate photon events from particle background. "Class 3" rejects the largest 
fraction of background contamination, and these events are used to study diffuse emission. We also require the zenith 
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angle be less than 105° to exclude most atmospheric gammas. By choosing to use these cuts, we study the events most 
likely to be real gamma rays, at the expense of a smaller effective area. These cuts discard roughly 3/4 of the signal, 
but vastly reduce the noise. 

The class 3 events are then binned in energy and into spatial pixels to produce counts maps. We use the hierarchical 
equal- area isolatitude pixelization (HEALPix), a convenient iso-latitude equal-area full-sky pixelization widely used in 
the CMB community.^'' 

Exposure maps 

Because the exposure on the sky is non-uniform, we generate an exposure map using the gtexpcube tool developed 
by the LAT team.^^ The exposure for each pixel is the LAT effective area at each 9 (angle with respect to the LAT 
axis) summed over the livetime of the LAT at that 9, and has units of cm^s. The exposure map spatially modulates 
the signal ±20% from raw photon counts and is slightly energy dependent for photon energies E > 1 GeV. To avoid 
systematic errors in generating the exposure maps we use the same energy grid used for the spectral plots, which is 
given in Table [TJ The default spatial binning of 1° in latitude and longitude is adequate, as neighboring bins have 
exposure differences of < 0.5%. The LAT contains 12 thin layers of tungsten designated "front" and 4 thicker layers 
designated "back." Photons may convert to an e+e" pair in any of the layers, and events are labeled "front" or 
"back" accordingly. The effective area as a function of energy is different for front and back events, so we use the 
"P6_V3_DIFFUSE : : FRONT" and "P6_V3_DIFFUSE : :BACK" instrumental response functions, respectively. Maps of counts 
are divided by the exposure and pixel solid angle to produce intensity maps (cts/s/cm^/sr). The "combined" maps 
are simply an exposure-weighted linear combination of the front and back intensity maps. The full-sky Fermi maps 
are displayed in Figure [T] along with the exposure map and mask. 

Smoothing 

In our analysis we wish to compare Fermi maps at different energies to each other, and to other templates. In order 
to match the PSFs of all these maps, we smooth each map by an appropriate kernel. On average, the front and back 
converting events have different PSFs (by a factor of ^ 2), so we smooth the intensity maps of each to a common PSF 
(usually 2°) to produce "front" and "back" smoothed maps (which can be averaged to obtain "combined" smoothed 
maps as above). 

For the PSF, we use a simple fit to the radius of 68% flux containment (in degrees): 

r68 = ((cii?^)"+c^)i/" (Bl) 

where E is in GeV, a = 1.2, /S = -0.83, and (ci,C2) = (0.50,0.04) and (0.90,0.09) for front- and back-converting 
events, respectively. This yields rgs = ciE^ at low E and rgs = C2 at high E. For simplicity, we assume the PSF is 
a Gaussian so that the FWHM, / = 2r5o ~ 1.56r68. The "raw" FWHM of a counts map is then taken to be Eq. IBll 
evaluated at i?mean = V^oEi, where Eq < E < Ei for the events in the map. In order to smooth to the target PSF 
(usually /targ — 2°), a Gaussian smoothing kernel of size /korn is used such that 

/kern = \J ftarg ~ /raw i^'^) 

For the lowest energy maps {E < 1 GeV), where /raw is large, we take /targ = 3 or 4°. 

Point Source Mask 

The point source mask (Figure [1} contains the 3-month lOcr point source catalog, plus the LMC, SMC, Orion, and 
Cen A. The mask radius for point sources is taken to be the 95% containment radius for the lowest energy event in 
the energy bin. For unsmoothed maps, both the counts map and the exposure map are multiplied by the mask. For 
smoothed maps, counts and exposure are multiplied by the mask, then smoothed, then divided. Pixels where the 
smoothed masked exposure drops below 25% of the mean unmasked exposure are replaced with zeros. In cases where 
the mask radius is much smaller than the smoothing radius, the mask is "smoothed away." In other cases, masked 
pixels are visible in the smoothed maps. Because of the energy dependence of the mask radius, the effective mask 
changes with energy, so care must be taken in cross comparisons between energy bins, especially within a few degrees 
of the Galactic plane. 

Data Release 

All LAT maps described in this section are available in the HEALPix pixelization on the web.^^ Maps of front- and 
back- converting events are available, as well as the combined maps. All three are available smoothed and unsmoothed, 
and with and without point source masking, for a total of 12 maps at each energy. Both the imbin and specbin energy 
binnings are available, for a total of 240 FITS files. Grayscale and color jpegs are also available to provide a quick 
overview. Software used to make the maps is available on request. 

HEALPix software and documentation can be found at http://healpix.jpl.nasa.gov, and the IDL routines used in this analysis are 
available as part of the IDLUTILS product at http://sdss3data.lbl.gov/software/idlutils. 
http: //fermi .gsfc .nasa.gov/ssc/data/analysis/documentation 
^^•^ http://feriiii.skymaps.info 
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Fig. 18. — Left panel: a 2° FWHM mock map fit to a 2° model. Right panel: a 1° mock map is smoothed to a 2° PSF and then fit 
to the 2° model. In each case, 1000 realizations are computed, and the histograms of x a-re shown for the four parameters. In each case, 
the X histogram has a mean near zero and RMS near 1, indicating that the error analysis is correct, even in the case of a smoothed map. 
Biases introduced are small (<g; Icr) for the case shown, and are relatively even less important at higher energy. 

POISSON LIKELIHOOD ANALYSIS 

Ideally the spatial templates (dust map, synchrotron map, etc.) going in to the Type 3 likelihood analysis (t ^3.3p 
would be smoothed to the PSF of the gamma-ray data, allowing an event-by-event likelihood to be calculated. However, 
at some energies, the LAT has a higher angular resolution than our templates. In this case, the LAT data and templates 
must be smoothed to a common PSF. In this section we explain the Poisson likelihood analysis in more detail, and 
check that a Poisson likelihood evaluated on a smoothed map behaves as expected. 

For the Type 3 fit we maximize the Poisson likelihood of the 4-template model in order to weight the Fermi data 
properly. In other words, for each set of model parameters, we compute the log likelihood 



ln£ = ^ [ki InyUj - i^i - In(fcJ)] , 



(CI) 



where /i^ is the synthetic counts map (i.e., linear combination of templates times exposure and mask) at pixel i, and k 
is the map of observed counts. The fj,i term depends only on the model parameters, and is not affected by smoothing 
the data. The ln(fci!) term does not depend on the model parameters, and so it cannot affect the best-fit model or 
uncertainties. Furthermore, there is no problem evaluating the likelihood for fractional k using the gamma function. 
So the ki In is the only potentially problematic term. 

In order to investigate the effects of smoothing, we generate each mock map by taking a linear combination of 
templates (including the uniform background) , multiplying it by the exposure and mask to obtain a map of predicted 
counts, and Poisson sampling it in HEALPix pixels. This map is then passed to our parameter estimation code to 
obtain the best fit values and uncertainties for the 4 parameters. Repeating this procedure for 1000 mock maps, we 
compute 

X=(^^) (C2) 

for each parameter and plot the histograms (Figure I18|) . An unbiased fit corresponds to (x) = and a correct 
uncertainty estimate corresponds to stdev(x) = 1- The uncertainties appear to be correctly estimated and the bias 
is small. Note that this behavior is different from an exp(— x^/2) likelihood on a smoothed map: in that case the 
correlated noise induced by the smoothing must be taken account of carefully. 

The above analysis was done for a range of energy bins with similar results. The analysis was also completely redone 
with a minimization instead (indeed, this provides the initial guess for the Poisson likelihood minimization) with 
very similar results. Because we feel the Poisson likelihood is more correct for the problem at hand, we use it for the 
"official" results. 



